Integrating evolution into ecological modelling: accommodating phenotypic changes in agent based models.

PMCID: PMC3733718This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.Evolutionary change is a characteristic of living organisms and forms one of the ways in which species adapt to changed conditions. However, most ecological models do not incorporate this ubiquitous phenomenon. We have developed a model that takes a 'phenotypic gambit' approach and focuses on changes in the frequency of phenotypes (which differ in timing of breeding and fecundity) within a population, using, as an example, seasonal breeding. Fitness per phenotype calculated as the individual's contribution to population growth on an annual basis coincide with the population dynamics per phenotype. Simplified model variants were explored to examine whether the complexity included in the model is justified. Outputs from the spatially implicit model underestimated the number of individuals across all phenotypes. When no phenotype transitions are included (i.e. offspring always inherit their parent's phenotype) numbers of all individuals are always underestimated. We conclude that by using a phenotypic gambit approach evolutionary dynamics can be incorporated into individual based models, and that all that is required is an understanding of the probability of offspring inheriting the parental phenotype

How predation and landscape fragmentation affect vole population dynamics

Background: Microtine species in Fennoscandia display a distinct north-south gradient from regular cycles to stable populations. The gradient has often been attributed to changes in the interactions between microtines and their predators. Although the spatial structure of the environment is known to influence predator-prey dynamics of a wide range of species, it has scarcely been considered in relation to the Fennoscandian gradient. Furthermore, the length of microtine breeding season also displays a north-south gradient. However, little consideration has been given to its role in shaping or generating population cycles. Because these factors covary along the gradient it is difficult to distinguish their effects experimentally in the field. The distinction is here attempted using realistic agent-based modelling. Methodology/Principal Findings: By using a spatially explicit computer simulation model based on behavioural and ecological data from the field vole (Microtus agrestis), we generated a number of repeated time series of vole densities whose mean population size and amplitude were measured. Subsequently, these time series were subjected to statistical autoregressive modelling, to investigate the effects on vole population dynamics of making predators more specialised, of altering the breeding season, and increasing the level of habitat fragmentation. We found that fragmentation as well as the presence of specialist predators are necessary for the occurrence of population cycles. Habitat fragmentation and predator assembly jointly determined cycle length and amplitude. Length of vole breeding season had little impact on the oscillations. Significance: There is good agreement between our results and the experimental work from Fennoscandia, but our results allow distinction of causation that is hard to unravel in field experiments. We hope our results will help understand the reasons for cycle gradients observed in other areas. Our results clearly demonstrate the importance of landscape fragmentation for population cycling and we recommend that the degree of fragmentation be more fully considered in future analyses of vole dynamics

A Software Tool to Model Genetic Regulatory Networks. Applications to the Modeling of Threshold Phenomena and of Spatial Patterning in Drosophila

We present a general methodology in order to build mathematical models of genetic regulatory networks. This approach is based on the mass action law and on the Jacob and Monod operon model. The mathematical models are built symbolically by the Mathematica software package GeneticNetworks. This package accepts as input the interaction graphs of the transcriptional activators and repressors of a biological process and, as output, gives the mathematical model in the form of a system of ordinary differential equations. All the relevant biological parameters are chosen automatically by the software. Within this framework, we show that concentration dependent threshold effects in biology emerge from the catalytic properties of genes and its associated conservation laws. We apply this methodology to the segment patterning in Drosophila early development and we calibrate the genetic transcriptional network responsible for the patterning of the gap gene proteins Hunchback and Knirps, along the antero-posterior axis of the Drosophila embryo. In this approach, the zygotically produced proteins Hunchback and Knirps do not diffuse along the antero-posterior axis of the embryo of Drosophila, developing a spatial pattern due to concentration dependent thresholds. This shows that patterning at the gap genes stage can be explained by the concentration gradients along the embryo of the transcriptional regulators

Host plant quality, spatial heterogeneity, and the stability of mite predator–prey dynamics

Population dynamics models suggest that both the over-all level of resource productivity and spatial variability in productivity can play important roles in community dynamics. Higher productivity environments are predicted to destabilize consumer–resource dynamics. Conversely, greater heterogeneity in resource productivity is expected to contribute to stability. Yet the importance of these two factors for the dynamics of arthropod communities has been largely overlooked. I manipulated nutrient availability for strawberry plants in a multi-patch experiment, and measured effects of overall plant quality and heterogeneity in plant quality on the stability of interactions between the phytophagous mite Tetranychus urticae and its predator Phytoseiulus persimilis. Plant size, leaf N content and T. urticae population growth increased monotonically with increasing soil nitrogen availability. This gradient in plant quality affected two correlates of mite population stability, population variability over time (i.e., coefficient of variation) and population persistence (i.e., proportion of plant patches colonized). However, the highest level of plant quality did not produce the least stable dynamics, which is inconsistent with the “paradox of enrichment”. Heterogeneity in plant productivity had modest effects on stability, with the only significant difference being less variable T. urticae densities in the heterogeneous compared to the corresponding homogeneous treatment. These results are generally congruent with metapopulation theory and other models for spatially segregated populations, which predict that stability should be governed largely by relative movement rates of predators and prey—rather than patch quality

Long-term cyclic persistence in an experimental predator–prey system

Predator–prey cycles rank among the most fundamental concepts in ecology, are predicted by the simplest ecological models and enable, theoretically, the indefinite persistence of predator and prey1,2,3,4. However, it remains an open question for how long cyclic dynamics can be self-sustained in real communities. Field observations have been restricted to a few cycle periods5,6,7,8 and experimental studies indicate that oscillations may be short-lived without external stabilizing factors9,10,11,12,13,14,15,16,17,18,19. Here we performed microcosm experiments with a planktonic predator–prey system and repeatedly observed oscillatory time series of unprecedented length that persisted for up to around 50 cycles or approximately 300 predator generations. The dominant type of dynamics was characterized by regular, coherent oscillations with a nearly constant predator–prey phase difference. Despite constant experimental conditions, we also observed shorter episodes of irregular, non-coherent oscillations without any significant phase relationship. However, the predator–prey system showed a strong tendency to return to the dominant dynamical regime with a defined phase relationship. A mathematical model suggests that stochasticity is probably responsible for the reversible shift from coherent to non-coherent oscillations, a notion that was supported by experiments with external forcing by pulsed nutrient supply. Our findings empirically demonstrate the potential for infinite persistence of predator and prey populations in a cyclic dynamic regime that shows resilience in the presence of stochastic events

Effects of Sample Size on Estimates of Population Growth Rates Calculated with Matrix Models

BACKGROUND: Matrix models are widely used to study the dynamics and demography of populations. An important but overlooked issue is how the number of individuals sampled influences estimates of the population growth rate (lambda) calculated with matrix models. Even unbiased estimates of vital rates do not ensure unbiased estimates of lambda-Jensen's Inequality implies that even when the estimates of the vital rates are accurate, small sample sizes lead to biased estimates of lambda due to increased sampling variance. We investigated if sampling variability and the distribution of sampling effort among size classes lead to biases in estimates of lambda. METHODOLOGY/PRINCIPAL FINDINGS: Using data from a long-term field study of plant demography, we simulated the effects of sampling variance by drawing vital rates and calculating lambda for increasingly larger populations drawn from a total population of 3842 plants. We then compared these estimates of lambda with those based on the entire population and calculated the resulting bias. Finally, we conducted a review of the literature to determine the sample sizes typically used when parameterizing matrix models used to study plant demography. CONCLUSIONS/SIGNIFICANCE: We found significant bias at small sample sizes when survival was low (survival = 0.5), and that sampling with a more-realistic inverse J-shaped population structure exacerbated this bias. However our simulations also demonstrate that these biases rapidly become negligible with increasing sample sizes or as survival increases. For many of the sample sizes used in demographic studies, matrix models are probably robust to the biases resulting from sampling variance of vital rates. However, this conclusion may depend on the structure of populations or the distribution of sampling effort in ways that are unexplored. We suggest more intensive sampling of populations when individual survival is low and greater sampling of stages with high elasticities

Optimizing Metapopulation Sustainability through a Checkerboard Strategy

The persistence of a spatially structured population is determined by the rate of dispersal among habitat patches. If the local dynamic at the subpopulation level is extinction-prone, the system viability is maximal at intermediate connectivity where recolonization is allowed, but full synchronization that enables correlated extinction is forbidden. Here we developed and used an algorithm for agent-based simulations in order to study the persistence of a stochastic metapopulation. The effect of noise is shown to be dramatic, and the dynamics of the spatial population differs substantially from the predictions of deterministic models. This has been validated for the stochastic versions of the logistic map, the Ricker map and the Nicholson-Bailey host-parasitoid system. To analyze the possibility of extinction, previous studies were focused on the attractiveness (Lyapunov exponent) of stable solutions and the structure of their basin of attraction (dependence on initial population size). Our results suggest that these features are of secondary importance in the presence of stochasticity. Instead, optimal sustainability is achieved when decoherence is maximal. Individual-based simulations of metapopulations of different sizes, dimensions and noise types, show that the system's lifetime peaks when it displays checkerboard spatial patterns. This conclusion is supported by the results of a recently published Drosophila experiment. The checkerboard strategy provides a technique for the manipulation of migration rates (e.g., by constructing corridors) in order to affect the persistence of a metapopulation. It may be used in order to minimize the risk of extinction of an endangered species, or to maximize the efficiency of an eradication campaign

Man Bites Mosquito: Understanding the Contribution of Human Movement to Vector-Borne Disease Dynamics

In metropolitan areas people travel frequently and extensively but often in highly structured commuting patterns. We investigate the role of this type of human movement in the epidemiology of vector-borne pathogens such as dengue. Analysis is based on a metapopulation model where mobile humans connect static mosquito subpopulations. We find that, due to frequency dependent biting, infection incidence in the human and mosquito populations is almost independent of the duration of contact. If the mosquito population is not uniformly distributed between patches the transmission potential of the pathogen at the metapopulation level, as summarized by the basic reproductive number, is determined by the size of the largest subpopulation and reduced by stronger connectivity. Global extinction of the pathogen is less likely when increased human movement enhances the rescue effect but, in contrast to classical theory, it is not minimized at an intermediate level of connectivity. We conclude that hubs and reservoirs of infection can be places people visit frequently but briefly and the relative importance of human and mosquito populations in maintaining the pathogen depends on the distribution of the mosquito population and the variability in human travel patterns. These results offer an insight in to the paradoxical observation of resurgent urban vector-borne disease despite increased investment in vector control and suggest that successful public health intervention may require a dual approach. Prospective studies can be used to identify areas with large mosquito populations that are also visited by a large fraction of the human population. Retrospective studies can be used to map recent movements of infected people, pinpointing the mosquito subpopulation from which they acquired the infection and others to which they may have transmitted it